Relaxation of a qubit measured by a driven Duffing oscillator 



O 

(N 



43 

^— > ■ 
c ■ 
ctf ■ 
=5 

cr 



(N 
> 

00 

o 
o 



X 



I. Serban, 1 ' 2 - 3 M. I. Dykman, 4 <E and F. K. Wilhelm 1 '0 

1 IQC and Dept. of Physics and Astronomy, University of Waterloo, 
200 University Ave W, Waterloo, ON, N2L 3G1, Canada 
1 Department Physik, Arnold- Sommerf eld- Center for Theoretical Physics, and Center for Nanoscience, 
Ludwig-Maximilians-Universitat, Theresienstr. 37, 80333 Munchen, Germany 
5 ' Instituut Lorentz, Leiden University, Niels Bohrweg 2, NL-2333 CA Leiden, The Netherlands 
^Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 

(Dated: January 11, 2010) 

We investigate the relaxation of a superconducting qubit for the case when its detector, the Joseph- 
son bifurcation amplifier, remains latched in one of its two (meta)stable states of forced vibrations. 
The qubit relaxation rates are different in different states. They can display strong dependence on 
the qubit frequency and resonant enhancement, which is due to quasienergy resonances. Coupling 
to the driven oscillator changes the effective temperature of the qubit. 
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I. INTRODUCTION 

The Duffing oscillator is a paradigmatic model for 
nonlinear dynamics. Its quantum properties [H-Q have 
gained renewed recent interest in the light of nanome- 
chanics [1] and as a qubit readout device [l0l - fl3j for 
superconducting qubits [Hj, the Josephson bifurcation 
amplifier (JBA). In the context of readout, examples of 
the experimental setup are a driven nonlinear supercon- 
ducting resonator [l3[ coupled to the quantronium qubit 
[lH or a superconductin g q uantum interference device 
(SQUID) used as a JBA [12|], inductively coupled to the 
flux qubit and in contact with dissipative measurement 
circuitry. 

The JBA detection profits from the very different am- 
plitudes and phases of the coexisting states of oscilla- 
tor's forced vibrations and the formally infinite gain at 
the bifurcation pjj| - [T3| . In the first step of detection, the 
JBA is driven through the bifurcation point where, de- 
pending on the state of the qubit, one of the attractors 
disappears. The JBA then decides between the attrac- 
tors thus pre- measuring the qubit. This decision is prob- 
abilistic, because of quantum and classical fluctuations 
that can cause an inter-attractor transition even before 
the bifurcation point is reached. At the next stage, the 
parameters are changed so as to suppress the probability 
of fluctuation-induced switching and the latched vibra- 
tion state is observed. Qubit energy relaxation in the 
latter state will drive the measurement away from the 
ideal quantum nondemolition regime in that the post- 
measurement state is not correlated with the measure- 
ment output [jq . 

In this paper, we study qubit relaxation mediated by 
the oscillator. Because of the oscillator-qubit interaction, 
coupling of the oscillator to a heat bath leads to an indi- 
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rect coupling of the qubit to that same bath and, conse- 
quently, to relaxation. The energy emitted by the qubit 
creates excitations in the bath. The qubit transition can 
be stimulated by forced vibrations of the oscillator, which 
have a comparatively large amplitude already for a weak 
resonant driving. This means that the driving accelerates 
qubit relaxation. For the lowest-order stimulated pro- 
cesses that we consider, the energy transfer to the bath 
becomes \u q ± ojf\, where io q is the qubit frequency and 
u>F is the driving frequency. As a consequence, the qubit 
relaxation rate becomes dependent on the vibration am- 
plitude of the oscillator, which bears on the experimental 
observations [Ta |. 



The dependence on the vibration amplitude is even 
more pronounced where u> q /2 is close to the oscillator 
eigenfrequency ujq, \uj q — 2wn| <C wo. In this case the 
decay rate displays resonances where h{uj q — 2wo) goes 
through the quasienergy level spacing of the driven oscil- 
lator. We have analyzed both resonant and nonresonant 
cases. The results coincide where the ranges of their ap- 
plicability overlap, which indicates that they describe all 
qubit frequencies. We have also extended the analysis to 
the case where uj q is close to ujq. 

In Sec. II we give a general expression for the relaxation 
rate of a qubit coupled to a driven quantum oscillator. 
In Sec. Ill we describe the theory of the oscillator and 
its fluctuations close to the stable vibrational states. In 
Sec. IV we provide an explicit expression for the oscil- 
lator power spectrum. The proposed method and the 
result go beyond the immediate problem of the qubit re- 
laxation. The obtained expression allows us to study the 
qubit relaxation rate as a function of frequency detuning 
uj q — 2bJF, including the onset of quasienergy resonances, 
and to discuss resonant heating and cooling of the qubit. 
In Sec. V qubit relaxation and the change of the effec- 
tive qubit temperature are studied far from resonance. 
Section VI provides a summary of the results. 
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II. GOLDEN RULE QUBIT RELAXATION 

We will express the qubit relaxation rate in terms of 
the power spectra of small-amplitude fluctuations of the 
driven oscillator about its steady states of forced vibra- 
tions. The fluctuations will be assumed small, which jus- 
tifies describing them in a linearized approximation. This 
approach allows us to find the relaxation rate both where 
the qubit frequency uj q is close or far detuned from the 
oscillator eigenfrequency ujq and 2wo- 

The Hamiltonian of a superconducting qubit reads 



K 5 ~ 



(1) 



where o~ x . y . z are Pauli matrices and \S\ <C w. The induc- 
tive coupling of the qubit to the oscillator will be treated 
as a weak perturbation and, for the typical setup of su- 
perconducting qubit measurements, is quadratic in the 
oscillator coordinate x [O, OJ, [H EH 



Hid = 7, a zX , 



(2) 



where A q is the oscillator frequency shift due to the cou- 
pling. We will also discuss the linear case, in which the 
operator x 2 in Eq. @ is replaced by x. As the interac- 
tion @ does not commute with the qubit Hamiltonian 
H q it introduces energy relaxation. 

We assume that the oscillator is coupled to the bath 
strongly enough, so that its damping k exceeds the decay 
rate of the qubit. Then the oscillator serves as a bath for 
the qubit. To second order in A q , the decay rate of the 
qubit excited state T e has a standard golden rule form 

r e = C r ReG«), C r = ^(mw A q S/hio q ) 2 , (3) 

G(u) = / dte lut (Sx 2 (t + t)Sx 2 {t)) f , 
Jo 



where u q = ^A" 2 



y/ w 2 + 5 2 is the qubit transition frequency 



and 5x 2 (t) = x 2 (t) — (x 2 (t)). The averaging (•) is done 
disregarding the qubit-oscillator interaction and assum- 
ing that the oscillator fluctuates about a given steady 
state of forced vibrations. In calculating (A(t + t)B{t))f 
we imply additional averaging with respect to r over the 
driving period; this period is small compared to T^ 1 due 
to the underdamped nature of the system. 

The expression for the excitation rate out of the 
qubit ground state T g has the same form as T e except 
G{uj q ) should be replaced by G(—uj q ). We will see that 
the ratio r ff /r e is not given by the Boltzmann factor 
exp(— Tuujq/kBT) due to the nonequilibrium character of 
this driven system. The Bloch-Redfield rates for the 
qubit are 



Tf 1 = r 

1 



•r„, 



(4) 



Of primary interest for us will be the decay rate Tj~ . 
As we show, it has a component where one quantum 
of the oscillator or the bath is created or annihilated; 
this component is proportional to the squared dimension- 
less amplitude r 2 [see Eq. [T0|) ] of the attractor the JBA 
has latched to; T^~ l also has a two-quantum component. 
Both of them can lead to enhanced relaxation at specific 
working points, identified as quasienergy resonances. 



III. DYNAMICS OF A RESONANTLY DRIVEN 
NONLINEAR OSCILLATOR 

To find the oscillator correlation function G(lu) we as- 
sume the limit where the resonant driving current Ib of 
the corresponding SQUID is much smaller than the criti- 
cal current I* but is large enough to reach the oscillator 
bistability regime. Then the driven SQUID can be mod- 
eled by a Duffing oscillator. Dissipation is described as 
resulting from contact with a bath of harmonic oscilla- 
tors. The Hamiltonian of the oscillator + bath reads 



H d {t) 
Hsif) 
Hp, 



H s {t) + H I + H Bl 



(5) 



p 2 musk „ 9 ^ d _ , . . 
-^x 2 - ls x A + F(t)x, 



2m 



H, 



1 2 



_ T -1 

2 Jl 



2C r (w/S) 2 ReG(0). 



where F(t) = F Q (e iUpt + e~ iujpt ) is the driving field with 
frequency up. The spectral density of the bath weighted 
with the interaction J(oj) = "X^j X 2 S(lu — ujj) is assumed 
to be a smooth function for ui close to ujq . The oscillator 
damping constant is k — J(wo)/2So;om (cf. [3]). For the 
case of an Ohmic bath we have J(w) = 2ftmKW0(w c — uj) 
where w c is a high-frequency cutoff. The interaction- 
induced renormalization of ujq is supposed to be incorpo- 
rated into wo- 

The steady states of the oscillator are forced vibra- 
tions. For weak damping and almost resonant driving, 
\ujp — luq\. k <C ojq, and up to moderate driving amplitude 
|JFb| these vibrations are almost sinusoidal. They can be 
analyzed in the rotating wave approximation (RWA) by 
changing to the rotating frame with a canonical transfor- 
mation U T (t) — exp(— iuptd^a) (a is the oscillator annihi- 
lation operator), and then by changing to dimensionless 
variables [!, 0, % 

x = C Ies (Q coscopt + Psmuipt), (6) 
p = — C res rrujjp(Q sinuipt — P coswpi), 

where C rcs = [2mwf (wp- ojf)/3js] 1 ^ 2 ■ This transforma- 
tion changes phase space volumes, so that 

[Q,P ]=iXs, A 5 = h Suj = uj f - w (7) 

2m z ijJ F |ocj| 

The dimensionless parameter Ag plays the role of 
Planck's constant in these units. We assume that As <C 
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1, in which case the oscillator dynamics is semiclassical. 
The oscillator has coexisting states of forced vibrations 
provided 7s (5w < 0. Motivated by the description of the 
JBA we have chosen the soft case of the Duffing oscilla- 
tor, 7s = towq/24 > 0, i.e., Su> < 0. The main results 
can be straightforwardly generalized to the hard case by 
choosing 5u > 0. 

If we disregard fast oscillating terms, following the 
RWA, the oscillator Hamiltonian Hs r = U]HsU r — 
ihU^dtUr becomes time-independent, 

H Sl = -(h\5u;\/\s)g, (8) 

where /3 is the dimensionless squared amplitude of the 
driving field, /3 = 3j s F$ /2(moj F \ 5lo \ f . 

For /? < 4/27 the Hamilton function Hs r has three ex- 
tremal points. In the presence of dissipation two of them, 
with largest and smallest Q, become attractors. In the 
laboratory frame they correspond to stable forced vibra- 
tions with amplitude C rcs r a . The third extremum of Hs r 
is a saddle point. For not too weak driving, the distance 
between the steady states in (P, Q) space is large com- 
pared to the typical amplitude of quantum fluctuations 

1 /2 

cx \g , as well as to the amplitude of classical fluctua- 
tions, which come into play for higher temperatures. 

Of primary interest for us are small-amplitude fluc- 
tuations of the oscillator about steady states. They 
can be conveniently described by the oscillator density 
matrix p in the Wigner representation, pw(Q,P) = 
J d£ exp[-i£P/A 5 ] {Q + i/2 \p\ Q - £/2) . When the oscil- 
lator is latched to one of the attractors, pw has a narrow 
Gaussian peak at that attractor. For weak damping, i.e., 
up to second order in the interaction Hj in Eq. ([5]), and 
linearized close to an attractor the master equation for 
p w in the RWA reads 0] 

p w w -V(K Pw ) +\ s k (n + V 2 Pw , (9) 

where we introduced vector K = 
(— \5w\dpg — kQ, \5uj\dQg — kP) that describes the 
drift of the mean coordinate; V = (<9q,<9p) [Hj], and 
h = l/[exp(fru)F/kE,T) — 1] is the Planck occupa- 
tion number of the oscillator. The attractor position 
(Qa,Pa) is found by demanding the drift to cease, 
K(Q a ,P a ) = 0, which leads to Q a = p- 1/2 r 2 a {r 2 a - 1), 
P a = -p-^K/lSuDfi, and 

rl[{rl-lf + {K/6uf] =P (10) 

(r 2 = Q 2 + P 2 ) J26[. Near attractor a vectors K should 
be expanded, 

K(Q, P) « [(ZV)K] a , Z = (Q - Q a , P - P a ). (11) 

Here, [•]„ indicates that the derivatives are calculated 
at attractor a. This expansion defines a matrix K. with 



elements IC nm — [d m K n ] a , where components 1 and 2 
correspond to Q and P, respectively. 

Equation (TlTj) is an expansion in the width of the Gaus- 
sian peak Q [As(2n+ l)] 1 / 2 as can be seen from combin- 
ing Eqs. ([9]) and (fTT]) . In our units, if this quantity is 
much smaller than unity, the deviations of Hs r , Eq. ([8} , 
from a quadratic form in P — P a and Q — Q a can be 
neglected. This is also reflected in the absence of higher 
derivatives and higher order terms in As in Eq. ©. 

IV. QUBIT DISSIPATION CLOSE TO 
RESONANCE WITH THE OSCILLATOR 

We first consider dissipation of the qubit in the sit- 
uation where its frequency is close to 2u>f sa 2ujq, i.e., 
\2cup — uj q \ however, we assume that \2cop — cu q \ 

largely exceeds the qubit dephasing rate, so that the 
driving does not pump the qubit resonantly. The ma- 
jor contribution to the correlation function G(uj q ) in 
Eq. Q comes from the terms in x{t\) proportional to 
[<9(ii) + iP(h)} exp(-iu F h) for ti = t + r and [Q(h) - 
iP(ti)]exp(iu F ti) for t\ — t (note that Q(t),P(t) re- 
main almost unchanged over a drive period). We can 
now re- write the excursion of the quadratic term as 

Sx 2 (t) = 2x a {t)[x(t) - x a {t)] + \x{t) - x a {t)} 2 . (12) 

Here, x a (t) is given by Eq. © with Q, P replaced by 
Q a ,P a , which have to be determined from Eq. (fT0|) . With 
x(t) — x a (t) being a Gauss distributed variable, Eqs. ^ 
and P^)l give rise to two types of decay rates. One is 
due to the noise of x(t) — x a (t) scaled by the prefactor 
x a . This is a one-quantum process, where a transition 
between the qubit states is accompanied by creation or 
annihilation of one quantum of the oscillator vibrations 
induced by the driving field, with the energy difference 
compensated by the field and the thermal bath. The 
second is based solely on the noise of [x(t) — x a (t)] 2 . This 
is a two-quantum process for the oscillator, which makes 
transitions between the next nearest levels with energy 
deficit compensated also by the field and the thermal 
bath. We are now going to compute both contributions. 

A. One-quantum noise 

We start with the first term in Eq. (TT2"]) , the one- 
quantum noise, i.e., the linear correlator scaled by the 
attractor radius. From Eqs. ^ 

GK)«^^r 2 iV + _K-2c, F ), (13) 

poo 

N+-(u)= / dtexp(iujt)(Z + (t)Z-(0)) alow , 
Jo 

where Z± = Z\ ± iZi and we define Z\ = Q — Q a , Zi = 
P — P a - The averaging (-^low is done in the rotating 
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frame and in the RWA. The time dependence of opera- 
tors is calculated with Hamiltonian Hs r + Hg + U^HiU r . 
This Hamiltonian has the same structure (except for the 
explicit form of Hg r ) as for a stationary oscillator with no 
driving, if one goes to the interaction representation via 
the transformation {7 r . Then we can apply the quantum 
regression theorem and write [24| , 



<Z„(i)Z m (0)> 



slow 



Tr[Z n p(t\Z„ 



(14) 



(n, m = 1, 2). The operator p(t\Z m ) is the weighted den- 
sity matrix; it satisfies the standard master equation [cf. 
Eq. J5J], but its initial condition is p(0\Z m ) = Z m p( st \ 
where p^ is the RWA density matrix in the steady state. 

Equation (flUf reduces the problem of the qubit decay 
rate to calculating the oscillator power spectrum in the 
rotating frame. For a quantum oscillator in a given 
attractor it can be done using Eqs. ©, (QUI) . (IT41 and 
the relations for the transition to the Wigner represen- 



ts 



tation Qp = (Q + ^i\ s d P ) pw{Q, P), Pp 
L J w L J 

(P- \i\ S d Q ) Pw {Q,P)- By multiplying Eq. © 

by Z n exp(iwi) and then integrating over Q,P,t 

we obtain a linear equation for matrix N(lu) of 

the Fourier-transformed correlators N nm (ui) = 

/ °° dtexp(iujt)(Z n (t)Z m (0)) s \ ovr , 



[iul + t)N{u)) = N - iiA s e, 
AA/C t + ftN = X s k (2n + 1) I, 



(15) 



where / is the unit matrix and e is the fully anti- 
symmetric tensor (the Levi-Civita tensor) of rank two. 
The matrix elements of the matrix Af are N nm = 
-j dQ dP Z n Z m p^ ] . From Eqs. (H3])-([T5]) one obtains 
for a given attractor a 



ReiV+_(w) = 2\ s n 



(n + l){[^-|H(2rg-l)] 2 + K 2 } + |H 2 
(w 2 - v 2 ) 2 + Ak 2 lo 2 

I 



(16) 



where v 2 a — n 2 + |<5w| 2 (3r^ — 4r 2 + 1) with r a obtained 
from solving Eq. 1101 In the high-temperature limit n ^> 1 
Eq. (fTo) coincides with the result obtained earlier [27[ in 
the classical limit. 



resonances for \uj q — 2u>f\ ~ 2v a , which will be larger 
than the background given by Eqs. (H2J), (fTo) for nr a < 
v a [Xs(2n+ 1)] 1 / 2 - Both their height and position depend 
on the field amplitude Fq. 



B. Two-quantum noise 



C. Discussion of the results for resonant decay 



The second contribution to the qubit decay rate orig- 
inates from direct two oscillator quanta transitions. It 
is not proportional to the squared attractor amplitude 
cx r 2 . The decay probability is determined by the 
power spectrum of the squared oscillator displacement 
[x(t) — x a (t)} 2 . At resonance, \w q — 2oj | <C uq, the 
main contribution to this power spectrum comes from 
the Fourier-transformed correlator (Z^(t)Z 2 (0)). The 
latter was studied earlier for a nonlinear oscillator in the 
absence of driving [24j and for linear driven oscillators 
[HjH^]. Up to moderate temperatures, hnjs /m 2 uj 2 K <C 
1, it leads to 



ReG(w g ) 



K (n+ l) 2 



(mwo) 2 (ui q — 2wq) 2 + 4k 2 



(17) 



This term becomes important for weak driving, where 
the amplitude of forced vibrations is on the order of the 
fluctuational oscillator displacement, r 2 < Xs(2n + 1). 
However, if the oscillator is strongly undcrdampcd in 
the rotating frame, k <C v a (see below), the decay rate 
may display attractor-dependent two-quanta quasienergy 



In total, the qubit decay rate is a sum of the one- 
quantum, Ti, and two-quantum rate, T2, i.e., r to t = Ti + 
T2- These contributions are described by Eq. ([3]) with 
the function G given by Eqs. (jl"3)) and (fTB]) for Ti and 
(HZD for r 2 . 

The one-quantum rate Ti is quadratic in the forced 
vibration amplitude cx r a . The physical origin of this 
term can be understood as follows: small fluctuations 
of x(t) — x a (t) are frequency-mixed with the forced vi- 
brations x a (t) in the laboratory frame, as described by 
the interaction cx Sx 2 . As a result fluctuation amplitude 
is multiplied by the forced-vibration amplitude, and the 
qubit decay rate is quadratic in the latter amplitude. 

Apart from the proportionality to r 2 , the attractor de- 
pendence of Ti is also due to the different curvature of 
the effective potentials around the attractors. For weak 
oscillator damping k <Si v ai the parameter v a in Eq. (|16[) 
is the frequency of small-amplitude vibrations about at- 
tractor a. It sets the spacing between the quasienergy 
levels, the eigenvalues of the rotating frame Hamiltonian 

Hsi close to the attractor. The function Re 7V_| (to) has 

sharp Lorentzian peaks at u) = ±v a with halfwidth k de- 
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tcrmined by the oscillator decay rate. The dependence 
of v a on the control parameter j3 is illustrated in Fig. [TJ 



1.2 




FIG. 1: Left panel: Squared scaled attractor radii r 2 a as func- 
tions of the dimensionless field intensity /3 for the dimension- 
less friction k/\5cj\ = 0.3. Right panel: The effective frequen- 
cies v a /\8ui\ for the same k/\Suj. Curves 1 and 2 refer to small- 
and large amplitude attractors. 

The decay rate of the excited state of the qubit T e oc 
Re (uj q — 2ujf) sharply increases if the qubit fre- 
quency ui q coincides with 2u F ± v a , i.e., ui q — 2lo f res- 
onates with the inter-quasienergy level transition fre- 
quency. This new frequency scale results from the in- 
terplay of the system nonlinearity and the driving and is 
attractor-specihc, as seen in Fig.[TJ In the experiment, for 
Lo q close to 2wo, the resonance can be achieved by tuning 
the driving frequency cu F and/or driving amplitude Fq. 
This quasienergy resonance destroys the QND character 
of the measurement by inducing fast relaxation. 

The analysis of the excitation rate out of the qubit 
ground state in the resonant case \u> q — 2luf\ -C oj f is sim- 
ilar; r g is given by Eqs. §5§ and (1151) with AT H — 2lo f ) 

replaced with N |_(— uj q + 2uj F ): the last function is given 

by Eq. (|16p with u — u> q — 2u F and with interchanged 
coefficients n + 1 O n. 

The scaled decay rates T e ^ g as functions of detuning 
uj q — 2u)p are illustrated in Fig. [5J Even for compara- 
tively strong damping, the spectra display well-resolved 
quasi-energy resonances, particularly in the case of the 
large-amplitude attractor. As the oscillator approaches 
bifurcation points where the corresponding attractor dis- 
appears, the frequencies v a become small (cf. Fig. [TJ 
and the peaks in the frequency dependence of r e ff move 
to uj q = 2u F and become very narrow, with width that 
scales as the square root of the distance to the bifurcation 
point. We note that the theory does not apply for very 
small u q — 2ujf\ where the qubit is resonantly pumped; 
the corresponding condition is (moJoA q 5C 2 cs rl/hujq) 2 <C 
[TiT 2 _1 + (uj q - 2uj f ) 2 T 1 T 2 }. For weak coupling to the 
qubit, T e -C k, it can be satisfied even at resonance. 

An important feature of the qubit relaxation in the 
presence of driving is that the stationary distribution 
over the qubit states differs from the thermal Boltzmann 
distribution. If the oscillator-mediated decay is the dom- 
inating qubit decay mechanism, the qubit distribution is 
determined by the ratio of the transition rates T e and 
T g . One can characterize it by effective temperature 
T e g = hoj q /[kB ln(r e /r g )]. If the term in curly brackets 
in the numerator of Eq. (|16[) is dominating, T e g w 2T, 




-3-2-10123 -3-2-10123 



FIG. 2: The scaled decay rate factors for the excited and 
ground states, curves 1 and 2, respectively, as functions of 
scaled difference between the qubit frequency and twice the 
modulation frequency; Fo = HCrra/6'ys- Left and right pan- 
els refer to the small- and large-amplitude attractors, with 
the values of /3 being 0.14 and 0.f2, respectively. Other pa- 
rameters are k/\5ui = 0.3, n = 0.5. 



but if the field parameters are varied so that this term 
becomes comparatively smaller T e g increases, diverges, 
and then becomes negative, approaching — 2T. Negative 
effective temperature corresponds to population inver- 
sion. The evolution of the effective temperature with the 
intensity of the modulating field is illustrated in Fig. [3J 




FIG. 3: The effective scaled qubit temperature T* s = 
kBT e s/hw q as function of the scaled field strength /3 in the 
region of bistability for the small- and large-amplitude at- 
tractors, left and right panels, respectively; {oJ q — 2wf)/|<5w| = 
—0.2 and 0.1 in the left and right panels; other parameters 
are the same as in Fig. [2] 



D. Linear qubit-to-oscillator coupling 

The above results can be extended also to a qubit lin- 
early coupled to the oscillator. Such coupling was dis- 
cussed in numerous contexts in the problem of a two-level 
system coupled to an oscillator; in particular, it underlies 
the broadly used Jaynes-Cummings model of the cavity 
quantum electrodynamics. The linear coupling leads to 
resonant decay of the qubit if the qubit frequency uj q is 
close to the oscillator eigenfrequency, \uj q — cjq\ <C uiq and 
the oscillator decay rate largely exceeds the qubit decay 
rate. The effect is significantly modified if the oscillator 
is resonantly driven. If we consider linear coupling of the 
Jaynes-Cummings form V x a x x, the corresponding contri- 
bution to the decay rate of the excited state of the qubit 
is 

T e = (V x w/hw q ) 2 (muj F \6u\/3 1 s)ReN + -(uq - lj f ) (18) 
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(for the coupling of the form V z a z x one should replace in 
the above expression V x w with V z 5). 

For linear coupling, the decay rate does not have the 
factor r 2 . However, it is still different for different at- 
tractors. If the oscillator is strongly underdamped, it 
displays resonances whenever h\u) q — up\ coincides with 
the oscillator quasienergy level spacing near the attractor 

hv a . 



V. QUBIT DECAY FAR FROM RESONANCE 

For the case where \u q — 2u$\ is not small compared 
to ujq and u q ^> v a , the correlation function G(oj q ) can 
be calculated by perturbation theory to the lowest order 
in the oscillator-to-bath coupling. To this end we again 
write Sx 2 (t) ss 2x a {t)[x(t) — x a (t)]. We need to find x(t) — 
x a {t) in the time range ~ \oj q — The analysis is 

familiar from the studies of phonon sidebands in solids 
[28j . We linearize the Heisenberg equation of motion for 
x(t) —x a (t). The terms oc 7s and oc F(t) in this equation 
can be disregarded, since the oscillator nonlinearity is 
weak, |7s|s 2 <C mu 2 , as is also the driving; the back- 
action (friction) force from the bath can be disregarded 
as well, since k <C \u q — 2uq\. However, the quantum 
force from the bath J2j + ~ b ]) 

may have components 
at frequency u q ±up and should be kept. Integrating the 
linear equation for x(t{) — x a (t\) first with t\ — t + r and 
then with t\ = r and substituting the result into Eq. ([3]), 
we obtain 



ReG(u q ) = 




where Ui takes on the values u q ± up, up — Lu q . The 
function J{uj) is the density of states of the bath weighted 
with the coupling to the oscillator, it is defined below 
Eq. and we assume J(u) = for u < 0. Function 
$r(wi) = n{uii) + 1 for uji = u q ±up and $x(u;i) = n{oJi) 
for Ui = up —u q \ here n(u) — 1 / [exp(?kj / ksT) — 1]. It is 
important that the two asymptotic expressions obtained 
in a different way, Eqs. (fTUf and (JTS]) on the one hand and 
Eq. (ITO)) on the other hand, coincide in the range v a , K <C 
\u q — 2up\ <C up, u c , thus indicating that we have found 
the qubit relaxation rates at arbitrary frequency. 

Equation (|I9[) describes decay of the qubit excitation 
into excitations of the bath coupled to the oscillator. This 
decay is mediated by the oscillator and stimulated by 
the driving. The decay rate is determined by the den- 
sity of states of the bath at the combination frequencies 
|w g ±wp|. Most notably, it is quadratic in the amplitude 
of the oscillator forced vibrations and therefore strongly 
depends on the occupied attractor. 

The far-from-resonance regime is important for the ex- 
periment (l6j with a SQUID-based bifurcation amplifier, 
since the amplifier was operated at a frequency far below 
the qubit frequency. In this experiment, the relaxation 
rate of the qubit for the large vibration amplitude of the 



oscillator was much larger than for the small amplitude, 
see Fig. 2c of Ref. [l|| and was increasing with the driving 
strength on the low-amplitude branch (branch 1 in the 
left panel of Fig. [I}, in qualitative agreement with the 
theory. It is not possible to make a direct quantitative 
comparison because of an uncertainty in the qubit relax- 
ation rates noted in Ref. [1(1 also, even though the decay 
rate of the oscillator at frequency u = uf is given, the os- 
cillator decay rate at the relevant much higher frequency 
u) = u q is not known. 

The excitation rate out of the ground state of the 
qubit T g is also determined by Eq. (1191) with interchanged 
n(ui) n(ui) in function If the oscillator- 

mediated decay is the dominating qubit decay mecha- 
nism, the effective qubit temperature T e g depends on the 
interrelation between the values of J(\u q ± up\). This 
goes beyond the temperature analysis in the resonant 
case because in the former, only J(up) entered into the 
relaxation rate whereas, for strong detuning, the environ- 
mental spectrum is probed at multiple, distinct frequen- 
cies. The analysis here is fully similar to that of cooling 
and heating of oscillato rs )24l.l29ll , the area that attracted 
much interest recently [30L 13 llj . We note that, in contrast 
to the previous work, here we consider heating and cool- 
ing of a two-level system, a qubit, that results from the 
coupling to an oscillator rather than heating and cooling 
of the oscillator that may result from such coupling pro- 
vided the relaxation rate of the qubit is higher than that 
of the oscillator. 

If J(u q — up) ^> J(u q + up), the qubit temper- 
ature T e g w Tu q /(u q — up) exceeds T, whereas for 
J(u q — up) J(u q +up) we have driving-induced cool- 
ing of the qubit, T e g « Tu q /(u q + up); for up > u q the 
temperature T c g w —Tu q /{up — u q ) is negative. Heating 
and cooling of the qubit can be achieved also for nonres- 
onant excitation of the oscillator [24, 29j, |30( • 

For weak driving and u q far from 2u$, qubit decay 
may be due to two-quanta transitions in which an oscil- 
lator makes a transition between its neighboring levels 
and a quantum of the bath is created or annihilated (for 
stronger oscillator-bath coupling, two bath quanta tran- 
sitions may be also important). The corresponding rate 
is determined by ReG(a; g ) « (2h/m 3 uo)J2 i J{ui)(uQ ~ 
uf ) _2 $T(wi)$Ti, where Ui — ±u q ± u . Here, func- 
tion $t(^) is the same as in Eq. (fT9")) (except that up 
should be replaced with uq), whereas $Ti = n + 1 for 
Ui = ± (u q — uq and <&Ti — n for Ui = u q + uq . This ex- 
pression coincides with Eq. (fl"7|) for uq ^ \u q — 2uq\ ^ k. 

The results on nonresonant oscillator-mediated decay 
can be extended also to the case of linear qubit-oscillator 
coupling, with coupling energy V x a x x. In this case 
r e = 2(V x w/hmu q ) 2 (u 2 q - u 2 ) l2 J{u q )[n(u q ) + 1]. The 
decay rate is independent of the attractor occupied by 
the oscillator. 

Due to the factor (w/S) 2 in T 2 _1 , i.e. the particular 
choice of working point, the dephasing rate of the qubit 
is several orders of magnitude larger that the relaxation 
rate, allowing for a good, near-QND, qubit detection. 



7 



VI. CONCLUSION 

In conclusion, the presented approach describes the re- 
laxation of a qubit coupled to an underdamped driven 
quantum oscillator, both in the regime where the qubit 
transition frequency u q resonates with the oscillator fre- 
quency or its second overtone and where it is far away 
from resonance. The proposed mechanism shows to what 
extent quantum measurements with bifurcation ampli- 
fiers are of the nondemolition type. 

The qubit relaxation rate is expressed in terms of the 
power spectrum of the oscillator. We find this power 
spectrum in the explicit form near resonance where, as 
we show, it can have a double-peak structure as function 
of frequency. This result and the proposed method are 
general for a nonlinear oscillator, and thus go beyond the 
problem of qubit relaxation. 

Both in the resonant and nonresonant regimes the 
qubit decay rate strongly depends on the attractor the 
oscillator is latched to. Far from resonance, including 
the experimentally important range of a high qubit fre- 



quency, and for the coupling quadratic in the oscillator 
coordinate the qubit decay rate is quadratic in the am- 
plitude of the oscillator's forced vibrations. For u! q close 
to 2loq, the decay rate is resonantly enhanced. For a 
strongly underdamped oscillator it displays narrow max- 
ima as a function of the control parameters once hu q goes 
through the appropriate oscillator quasienergy, which de- 
pends on the attractor. We observe that the decay me- 
diated by a driven oscillator changes the qubit temper- 
ature, leading to heating, cooling, or population inver- 
sion depending on the frequency detuning. In the case 
of linear qubit-oscillator coupling, the qubit relaxation 
rate does not contain the squared vibration amplitude as 
a multiplying factor, but for u q close to ujq it displays 
attractor-dependent quasienergy resonances. 
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